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A braneworld model for neutrino Dark Energy (DE) is presented. We consider a five dimensional 
two-branes set up with a bulk scalar field motivated by supergravity. Its low-energy efi'ective theory is 
\^ • derived with a moduli space approximation (MSA). The position of the two branes are parameterized 

' by two scalar degrees of freedom (moduli). After detuning the brane tensions a classical potential 

, for the moduli is generated. This potential is unstable for dSi branes and we suggest to consider as 

. a stabilizing contribution the Casimir energy of bulk fields. In particular we add a massive spinor 

(neutrino) field in the bulk and then evaluate the Casimir contribution of the bulk neutrino with the 
I help of zeta function regularization techniques. We construct an explicit form of the 4D neutrino 

mass as function of the two moduli. To recover the correct DE scale for the moduli potential 
the usual cosmological constant fine-tuning is necessary, but, once accepted, this model suggests a 
^>f^ I stronger connection between DE and neutrino physics. 
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I. INTRODUCTION 



There is an increasing observational evidence that our Universe is in a state of accelerated expansion fl] . The mixing 
' of data coming from CMB, supernovae la and abundance measurements illustrate a surprising picture of the cosmic 
energy budget: the dominant energy component (dark energy, DE, for reviews see [2|) has negative enough pressure 
and its contribution is comparable, at present, with that of the pressureless dark mSter component (DM), although 

■ their dynamical evolution during the cosmological history can be totally different ("coincidence problem"). To remove, 
, at least partially, this problem, a direct interaction between DM and DE has been discussed [3, 4, 5, 6, 7] and the 

IL" ' "old" idea of varying-mass particles [81 has been taken into account again. These models are alsPwSlSSi'rated from 
. 5^ I theoretical point of view by the identftncation of the DE with the string-theory dilaton f [9. 10^.^1 IJ . but also [12J for DE 
. as a stabilized dilaton). These so-called "coupled DE models" can show a late-time attracxor solucion, i.e. the late-time 
\^ ' cosmology is insensitive to the initial conditions for DE. For a phenomenologically acceptable scenario, in these models 

■ the present coupling of the scalar field with baryons must be sufficiently weak while the couplings with DM must be 
sufficiently strong (baryo-phobic quintessence) . This is a very general problem: allowing a direct interaction between 
matter and a quintessence (ultra-light ~ 10~'^^eV) field could be phenomenologically dangerous (violations of the 
equivalence principle, time dependence of couplings - for reviews see [13]'). 

A possible way-out of these problems is guaranteed by scalar-tenSnr (ST) theories of gravitation [14|, where by 
construction matter has a purely metric coupling with gravity (for reviews see [15] . [16] .[17J'). The existence of an 
attraction mechanism of ST towards Einstein's General Relativity (GR) has been atecuSSecTm [18. 19. 20] . Simultane- 
ously the presence of a second attractor is possible ([20]) and the correct evolution of poE aloitgci^-CcCned 'tracking' 



[21] solution is rescued. However deviations from standard cosmology and General Relativity may be present at nucle- 
osynthesis, CMB and solar system tests of gravity [22, 23]. Another possible way-out is to consider "chameleon fields" 
([24|): the mass of the scalar fields is determined bytnfeTiensity of the environment. On cosmological distances, where 
tfT?densities are very tiny, the mass can be of the order of the Hubble constant and the fields can roll on cosmological 
time scales. On the Earth instead, the density is much higher and the field acquires a large enough mass to evade 
experimental bounds mantaining relevant coupling with matter. 

Recently Fardon, Nelson and Wcincr (FNW) exploited the notion of "tracking" in a completely new way ([25]): 
they suggested that DE tracks neutrinos. They established a (not so!) firm connection between the dark energy 'scale 
and the neutrino mass scale: the similarity between the two scales becomes a true relation. The authors imagined 
furthermore that DE and neutrinos are directly coupled: a new dark sector arises, taking contributions from a scalar 
field (f) (a sort of quintessence) and neutrinos. In this scenario the dark energy density depends on the neutrino masses, 
which are not fixed but varying as functions of the neutrino density. The coupling between DE and neutrinos is encoded 
in the variable neutrino mass 'm^{(j)). During the cosmological expansion the temperature T decreases. When T < nii, 
neutrinos become non-relativistic and their contribution to the energy density becomes mass-dependent. In this way 
the coupling cj) — v brings to an effective potential of the form: 



K// = Vo{(j)) + n^m^(0), (1.1) 
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where Vq is the scalar field potential and is the neutrino number density. If we assume that Vq is minimized for 
large m^, a competition between the terms in is created, a minimum in Ve// arises and the neutrino mass is 

determined dynamically during the cosmological expansion. For works on mass varying neutrinos see also [26] , [27] . 
Several points are unsatisfactory about this approach (see also 

• The FNW model is characterized by Vo{(f)) (the potential for the scalar field) and m(0) (the neutrino mass). 
The explicit form of these functions is inserted by-hand. It would be interesting to obtain a firmer theoretical 
derivation of these functions in a "top-down" scenario. 

• The connection between dark energy and neutrino mass is not well established unless the scalar field gets an 
appropriate value. 

As we will show, an extradimensional scenario addresses both drawbacks of the FNW model: 

1) a 4D neutrino mass as a function of two scalar degrees of freedom (moduli) and the potential for these scalar fields 
are not inserted by hand but rather explicitly calculated. 

2) A firmer relationship is established between neutrino mass scale and Dark Energy scale. 

It must be stressed that solving these problems led us to a model that is significantly different: the only dynamic 
that we contemplate for the scalar fields is a chameleontic behaviour for one of the two moduli (the radion). Then, 
we have to distinguish between different distance scales: on cosmological distances, the densities are tiny and the 
radion can roll on cosmological time scales (the neutrino mass is variable with time and possible deviations from a 
pure cosmological constant contribution can be observed in the DE sector); on smaller scales instead, the chameleon 
acquires a large mass (radion stabilization) and both moduli are simply sitting on the minimum of the potential (for 
a discussion of possible manifestations of DE at short scales see [2^)- As a final result, the connection between DE 
and neutrino mass will be established in the form ' — ' 

( )mm ~ 1, (1.2) 

m„ 

where V is the total moduli potential. 

Before getting into the detailed presentation of the model, here is a summary of our approach. We consider a 
supersymmetric extension of Randall-Sundrum two branes model {\29, 30], while for reviews on extra-dimensions see 
[31]). This kind of approach finds its theoretical firmer justification in siring theory (see [32. 33j ) . Our general set-up 
wffl follow [34. 35] : a five-dimensional braneworld model with two boundary branes aiirna^calar field in the bulk 
motivated iVstrpergravity. As we will see in the forecoming sections, the most distinguishing feature of this (BPS) 
class of models is that the brane and bulk dynamics are related to each other and under control. Once we have solved 
bulk equations (Einstein equations and Klein-Gordon equation) the boundary conditions do not give any additional 
information, so we can put the branes without any obstacle into the background. The positions of the two branes 
are specified by two independent constants (no-force condition; for a discussion on BPS condition and connections 
between brane physics and black- hole physics see [36]). We will discuss than the moduli space approximation to 
derive the low energy action and a bi-scalar-tensor Ttleory is recovered, with the two scalar fields A and (moduli) 
parametrizing the positions of the branes: the two constants are promoted to scalar degrees of freedom. To render 
this model more realistic, we can leave the BPS configuration in two different ways: i) putting matter on the branes; 
ii) perturbing the brane tensions (i.e. parametrizing the supersymmetry breaking in a phenomenological way). It 
has been shown ([371) that it is possible to stabilize both moduli putting matter on branes without perturbing the 
brane tensions. HCTWever, since we want to describe a Universe with a positive cosmological constant, the detuning of 
the brane tensions is unavoidable and a classical potential for the moduli is generated ([38. 39] ). Unfortunately the 
potential is unstable and, to stabilize the system, we suggest to consider as a stabilizatiorr mechanism the quantum 
(Casimir) contribution from bulk fields. We will then add a massive spinor (neutrino) in the bulk and we will evaluate 
the Casimir energy using zeta function regularization techniques (for reviews on Casimir energy see [40] . [42] : for 
zeta function techniques see [41j, [42J, [43J). Remarkably the potential V ~ V{\,(t>) is not inserted by^tl&ncTtn our 
model and the geometrical colinguration of the system is stabilized (a long-standing problem in extra-dimensional 
and string-inspired scenarios, see for example [44] and references therein for a general treatment, but also [37, 45] for 
moduli stabilization in BPS braneworlds). — — 

As far as the 4d neutrino mass is concerned, following [46], we will obtain a small 4d neutrino mass without recurring 

to a see-saw mechanisrn. ... . 

About the organization of this article, in section 2 we present the action and the moduli space approximation; in 
section 3 we study the chameleon mechanism for the radion; in section 4 the explicit form of the neutrino mass is 
derived as function of the moduli; in section 5 the moduli potential is analyzed. In section 6 we will discuss the 
parameters of the model and the moduli stabilization. In section 7 we will draw some concluding remarks. The 
appendix contains some technical details. 



II. THE MODEL 



The bulk action is the sum of two main contributions: the first one (SUGRA) is motivated by supergravity and 
follows the second one (SPINOR) is the bulk- neutrino part. We have 

•Sbulk = 5'sUGRA + •S'sPINOR (2.3) 
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where the SUGRA action consists of two terms which describe gravity and the bulk scalar field (C) dynamics: 

1 /• , /„ 3 



5suGRA = 2^ y c^'^V^ (n-^ {{dCf + U)j; (2.4) 

while the SPINOR action, which takes into account a Dirac bulk fermion with mass m of order the fundamental 
scale M, is written in the form: 

S = Jd^Jdz^.f^E^ '-^j-{dA~^)'i'+'^^{r,a''}^ ~msgn(«)§vl/| . (2.5) 

We use capital indices A,B, . . . for objects defined in curved space, and lower-case indices a,b, . . . for objects defined 
in the tangent frame. The matrices 7° = (7^,i75) provide a four-dimensional representation of the Dirac matrices in 
five-dimensional flat space. The quantity is the inverse vielbein and wtcA is the spin connection. Because in our 
case the metric is diagonal the connection gives no contribution to the action (2.5). 

The sign change of the mass term under a — > —a is necessary in order to conserve a-parity, as required by the Z2 
orbifold symmetry we impose. 

Further, our setup contains two branes. One of these branes has a positive tension, the other brane has a negative 
tension. They are described by the action 



•S'branol = ~ TT^ / d X^/^^UbS{zi) , (2.6) 
5'brane2 = [ - QbU b6{z2) ■ (2.7) 

In these expressions, zi and Z2 are the (arbitrary) positions of the two branes, Ub is the superpotential; U ^ the bulk 
potential energy of the scalar field, is given by (BPS relation) 



This specific configuration, when no fields other than the bulk scalar field are present and when the bulk and brane 
potentials are unperturbed (no susy breaking effects), is the BPS configuration ([371). When Ub is the constant 
potential, the Randall- Sundrum model is recovered with a bulk cosmological constant A5 — (3/8) [/ = — (3/8)[/^. 
We will also include the Gibbons-Hawking boundary term for each brane, which have the form 

Sgh^ ^ f d'^xV^K, (2.9) 



where K is the extrinsic curvature of the individual branes. 

The solution of the system above can be derived from BPS-like equations of the form 



^__Ub ^,_dUB , . 



where ' = d/dz for a metric of the form 



ds^ = dz^ + a^{z)r]^^dx''dx''. (2.11) 
We will particularly focus on the case where the superpotential is an exponential function: 

Ub = 4fce"^. (2.12) 



The solution for the scale factor reads 



while the scalar field solution is 



In the a — > we retrieve the AdS profile 



a(z) = (1 - 4Wz)3^, (2.13) 



C ^ --\n(l~ Aka^z) . (2.14) 
a 



a(z)=e-'=^ (2.15) 



In that case the scalar field decouples and the singular point (for which the scale factor vanishes) at z» = l/(4fca^) is 
removed (singularities in braneworld scenarios are analyzed in f47]). 
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In the following we will discuss the moduli space approximation. Two of the moduli of the system are the brane 
positions. That is, in the solution above the brane positions are arbitrary. In the moduli space approximation, these 
moduli are assumed to be space-time dependent since in connection to matter distribution on branes. We denote the 
position of ("visible") brane 1 with zi = 4>{x'^) and the position of ("hidden") brane 2 with Z2 = A(a;^). We consider 
the case where the evolution of the brane is slow. This means that in constructing the effective four-dimensional 
theory we neglect terms like {d(f>)^. 

In addition to the brane positions, we need to include the graviton zero mode, which can be done by replacing 77^^ 
with a space-time dependent tensor gp_i,[x^). Note that the moduli space approximation is only a good approxima- 
tion if the time-variation of the moduli fields is small. This should be the case for late-time cosmology well after 
nucleosynthesis, which we are interested in. 

A. Moduli Space Approximation: the gravitational sector 

Let us first consider the SUGRA action. Replacing -q^i, with g^i,j{x^^) in (|2.11|l we have for the Ricci scalar TZ = 
IZ'^^^a? + TZ, where TZ is the Ricci-scalar of the background (|2.11|l . We explicitly use the background solution (|2.13|) 
and 12.14|l . so that R will not contribute to the low-energy effective action. Also, in this coordinate system, where 
the branes move, there is no contribution from the part of the bulk scalar field. Collecting everything we therefore 
have 



Shn 



Ik 



24 



with 



/(0,A) 



24 



dza?{z). 



(2.16) 



(2.17) 



We remind the reader that a{z) is given by (|2.13|l . 

We now turn to the boundary terms. The integrals H2.6|l and H2.7|l do not contribute to the effective action for 

the same reason that R does not contribute. Let us therefore turn our attention to the Gibbons-Hawking boundary 
terrns. 

First, it is not difficult to construct the normal vectors to the brane: 



1 



v/1 + (90)V«' 
Then the induced metric on each brane is given by 

inc 



1) • 



9 LiU 



Qd,2 



Thus, 



■9 



ind,l _ „4/ 



So the Gibbons-Hawking boundary terms take the form 



1 



2a2(0) 
1 

2a2(A) 



{dXf 



(2.18) 



(2.19) 
(2.20) 



(2.21) 
(2.22) 




1 



2o?{(j)) 
1 

The trace of the extrinsic curvature tensor can be calculated from 

1 

K -- 



{dXf 



K, 
K. 



(2.23) 
(2.24) 

(2.25) 



Neglecting higher order terms this gives 



1 - 



4a2 



(2.26) 
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The terms for the second brane can be obtained analogously. Using the BPS conditions and keeping only the kinetic 
terms, we get for the Gibbons-Hawking boundary terms 



+ 



3 

3 
4^ 



d'^xy/^a^ {X)Ub (A) (9(t)2 . 



Collecting all terms we find 
Smsa - 



'94. 



4t Kc: 4t hit: 



(2.27) 
(2.28) 

(2.29) 



Note that the kinetic term of the field has the wrong sign. This is an artifact of the frame we use here. It is possible 
to go to the Einstein frame with a simple conformal transformation, in which the sign in front of the kinetic term is 
correct for both fields. We observe that for a BPS system the moduli potential is totally absent: the fields describe 
flat directions. 

Coming back to the action above, we redefine the fields in the following way: 



2/3 



with 



then, the gravitational sector can be written as 



5'msa 



2fcK2(2a2 



1 



d'^x. 



/9 = 



-.94 



2a' 



4a2 



(2.30) 
(2.31) 

(2.32) 
(2.33) 



This is an action of the form of a multi-scalar tensor theory, in which one scalar field has the wrong sign in front 
of the kinetic term. Furthermore, in this frame there is a peculiar point where the factor in front of TZ can vanish, 
namely when <^ = A, which corresponds to colliding branes. 

In order to avoid mixed terms like {d^(j)){d^X), we shall define two new fields [831: 



(p = Qcoshi?, 
X — Q sinh R. 

To go to the Einstein frame we perform a conformal transformation: 

9tiu — Q 9tJ.v 

Then 



-9Q''n 



.9(7^-^(aQ) 



Collecting everything we get the action in the Einstein frame (where we now drop the tilde): 



-Sep — 



1 



2fcK2(2a2 



d^Xy 



-9 



1 + 2a' Q2 
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2a' + 1 



{dRf 



(2.34) 
(2.35) 



(2.36) 



(2.37) 



(2.38) 



Clearly, in this frame both fields have the correct sign in front of the kinetic terms. Note that for a ^ (i.e. the 
Randall-Sundrum case) the Q-field decouples. In this case, the field R plays the role of the radion, i.e. it measures 
the distance between the branes. Furthermore, we can identify the gravitational constant: 



IQttG = 2kKl{l + 2a'). 



(2.39) 
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B. Moduli space approximation: the matter sector 



In this section we are going to leave the BPS configuration. Let's suppose we put some matter on the branes: 
Standard Model fields on one brane ("visible") and Dark Matter particles on the other brane ("hidden"). When the 
cosmological evolution of the branes is analyzed in the absence of supersymmetry breaking potentials, it has been 
shown (\37]) that the branes are driven by their matter content: one brane is driven towards the minimum of Ub, 
while thfe~Other brane towards the maximum. In this way it is possible to achieve a full stabilization of the system, 
granted that we choose the superpotential of the theory in a proper way. If we don't localize matter on branes and 
don't break supersymmetry, than the BPS configuration is respected, the no-force condition between the branes is 
present and the system is static. 

However it is our intention to describe a Universe with a positive cosmological constant and detuning is necessary 
in order to achieve de Sitter geometry. We therefore introduce matter as well as supersymmetry breaking potentials 
V{Q, R) and W{Q, R) on each branes. We begin with the potentials: to first order in the moduli space approximation 
we get 

" d^xV^4[a\^)V{4>)] (2.40) 

with a-^icf)) = ,^4/(1+2"'). The expression for a potential W on the second brane is similar with a(0) is replaced by 
a(A). In the Einstein frame we have (dropping the tilde from the metric): 

rf^a;V^Q-»"'/(i+2"')(coshi?)4/(i+'"'V(g,i?) = ( d^x^V,s{Q,R), (2.41) 



where we have defined 

K=ff (Q, R) = g~^"'/('+'"')(coshi?)4/(i+2«')T/(Q, R). (2.42) 
The expression for W{Q, R) in the Einstein frame is 

j d^x^Q-'"'"'^^+^°'"\&mhRf/^^+^''"^W{Q.,R)= J d^x^W,s{Q, R), (2.43) 

where 

Weff(Q,i?) = g-«"'/(i+2"')(sinhi?)4/(i+2"')VK(Q,i?). (2.44) 
For matter, the action has the form 

=5W(vI/i,5»<^'i) and S^^^ ^ S^^ {^2, 9^11''% (2.45) 

where (7™'' denotes the induced metric on each branes and 'i'i the matter fields on each branes. Note that we do not 
couple the matter fields to the bulk scalar field, and thus not to the fields Q and R. In going to the Einstein frame 
we get 

S^r^^ = Si]^ {^i,A'{Q, R)g^,) and S^^ = S^^ {^>2, B\Q, R)g,,,), (2.46) 

In this expression we have used the fact that, in going to the Einstein frame, the induced metrics on each branes 
transform with a different conformal factor, which we denote with A and B. We have neglected the derivative terms 
in the moduli fields when considering the coupling to matter on the brane. They lead to higher order operators which 
can be easily incorporated. The energy-momentum tensor in the Einstein frame is defined as 

Ti^J^,^^S^rk\^fiQ,R)g) 

with an analogous definition for the energy-momentum tensor for matter on the second brane. 
In the Einstein frame, the total action is therefore, 

I2a^ {dQf 6 , 
l + 2a2 g2 2a2 + l^ ^ 



^'^^ = 16^ ' ^ 



- j d^x^iVM.R) + W,n{Q,R)) + S^^\-^^,A\Q,R)g^,) + S^'^\-^2.B''{Q,R)g^,). (2.48) 

We point out that in this braneworld scenario we naturally obtain a multimetric theory. Particles localized in 
different branes "feel" the respective induced metric {g^^^)- In the next section we will discuss an interesting phe- 
nomenological application of this fact: the chameleon mechanism. The moduli potential will be presented in a 
forecoming section. For more details on the moduli space approximation see [37] . 
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III. CHAMELEON FIELDS 



As we already mentioned in the introduction, the couphng of the moduU with matter is severely constrained by 
phenomenological tests. Introducing massless or ultra-light scalar fields with generic coupling with matter is prohibited 
by tests on the equivalence principle and on variations of the fundamental constants. Different ways to avoid this 
conflict have been considered in literature. Basically the fields introduced in the model have cither a small coupling 
today or they have a large enough mass (stabilization) . Recently a new stabilization scenario has been introduced in 
[24] . According to these authors, the mass of the scalar fields is determined by the density of the environment. On 
CoSnological distances, where the densities are very tiny, the mass can be of the order of the Hubble constant and the 
fields can roll on cosmological time scales. On the Earth instead, the density is much higher and the field acquires a 
large enough mass to evade all current experimental bounds on deviations from GR. In other words, the scalar field 
shows a " chameleontic behaviour". In this section we will discuss this kind of scalar fields. We will start by reviewing 
the mechanism with a basic example; then we will discuss the role of the radion as possible chameleon field. 



A. The chameleon mechanism 

Let's consider a multi-metric model with the following scalar-tensor action: 

S^J (3.49) 

where is a generic scalar field with potential V{(l)), -tjjm denotes matter fields which follow geodesies of a metric 
g^Iu = A^{<t')9tJ.v- The scalar field interacts with matter through the conformal factor ((/)), so the dynamics of the 
chameleon cj) is described by an effective potential that for pressureless (non-relativistic) matter takes the form (with 
general multi-metric theory a sum over i is present) 

V,s{4>)^V{4>)+ p^A{cp). (3.50) 

In other words, the conformal coupling modifies the Klein-Gordon equation for the chameleon in the following way: 

V^^ = a0p™A(0) , (3.51) 

where pm is conserved with respect to the Einstein frame metric (.g^jy). Supposing a run-away behaviour for the 
"bare" potential 1^(0), the effective potential will show a minimum if a "competition" between V{4)) and the density- 
dependent term is realized (i.e. A{(t)) is an increasing function). In this way the chameleontic properties of the field 
are recovered: the location of the minimum and the mass of the fluctuactions iv? — V^'j^ both depend on pm- 

Another effect which suppresses the chameleon-mediated force is the following: for large enough objects, the (/)-force 
on a test particle is almost totally due to a thin shell just below the surface of the object, while the contribution of 
core-matter is negligible. Only a small fraction of the total mass affects the motion of a test particle outside. This 
mechanism has been named "thin-shell mechanism". To give a more detailed derivation, we consider a chameleon 
solution for a spherically symmetric object of radius R and homogeneous density p. We take V((j)) = iW^+^/c/)", where 
M has units of mass, and an exponential coupling of the form A{(f>) = e^'^ with P — 0(1). Next we choose boundary 
conditions requiring that the solution be non-singular at the origin and that the chameleon tends to its environment 
value 00 far from the object. 

For sufficiently large objects the field assumes inside the object a value (j)c which minimizes the effective potential, 
i.e. V^{(f>c) + PpcC^'^"/^'^'''^ /Mpi = 0. This holds everywhere inside the object except within a thin shell of thickness 
Ai? below the surface where the field grows. Outside the object, the profile for </) is essentially that of a massive 
scalar, ~ exp(— mor)/r, where mo is the mass of the chameleon in the ambient medium. 

The thickness of the shell is related to 4>o, and the Newtonian potential of the object, $jv = M/SirMpiR, by 



Ai? 



R epMpi^N 
The exterior solution can then be written explicitly as [24] 



(3.52) 



- - (ii^ j [-JT ) ^ + ^0 . (3.53) 

This equation gives the correction to Newton's law at short distances in the form F ~ {1 + 9)Fm where = 
20^ (2Aii), Hence a thin-shell AR/R « 1 guarantees a small deviation from Newton's law. 
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B. The radion as a chameleon 



In the general set-up we described in the first section, the possible interpretation of the radion R (interbrane 
distance) as a dark energy chameleon field has been investigated in [481 . These authors introduced by-hand a run- 
away bare potential for the radion, to avoid the falling of the hidden tTrane into the naked singularity, in the form 

ViR) = A^R--'. (3.54) 

Their investigation studied the efficiency of the thin-shell mechanism for the radion. The coupling to matter was 
written by the expansion 

MR) ^1 + 1^ + -, (3.55) 

where the neglection of higher-order terms was allowed by the cosmological evolution, that drove the radion to small 
field value. Furthermore, the authors calculated a radion-mediated force, with a correction to Newton's law given by 

0=^Rl- (3.56) 

The strength of the radion-mediated force is specified by the cosmological value of the radion. Since 9 must be 
small the parameters of the theory must be adjusted to guarantee a small field value. In particular, imposing the 
gravitational constraints implied A ~ 10^"^ eV and 7 < 10~^. In other words the model must be extremely fine-tuned 
if we want simultaneously (1) to interpret the radion as a dark energy candidate and (2) to satisfy local gravitational 

constraints. The fine-tuning of 7 can be avoided by considering the potential V{R) — A''e^«^", but the fine-tuning 
on A remains as a cosmological constant fine-tuning. 

We are going to exploit the chameleon mechanism to stabilize the radion (R-modulus) following the approach of 
[481 . A detailed description of the moduli stabilization will be given in section 6. 

— To proceed further, our first step will be the creation of a "naturally small" neutrino mass in this braneworld set-up 
(see next section). The second step will be the construction of the potential for the scalar field(s) using a physical 
mechanism and not simply including potentials by- hand (see section 5). 



IV. THE NEUTRINO MASS 



In this section we describe how the sterile bulk fermion can provide a mechanism (different from the see-saw) to 
obtain a small Dirac neutrino mass in four dimensions. In our scenario all standard model matter and gauge fields 
are localized on the visible brane. We begin with a brief review of the scenario with mass- varying neutrinos recently 
proposed in [251. We proceed with the description of the Kaluza-Klein reduction, then we show that it is possible to 
localize the sTErile fermion near the hidden brane following [46] . In the end of the section we will discuss the role of 
the conformal transformation to the Einstein frame. 



A. Mass varying neutrinos and dark energy: the need for moduli stabilization 

We are now going to describe in more detail the idea proposed recently by Fardon, Nelson and Weiner [25] that 
connects the dark energy sector with neutrino's one. ' 

The main assumption in the FNW scenario is that the dark energy sector has two components: the dark energy 
contribution of some (for example) quintessence field (j) plus the neutrino energy density 

Pdark Pi^ + Pdark energy. (4.57) 

An interaction is present between the two components, due to the dependence of the neutrino mass on the 
quintessence field. The form of this dependence is obtained via a see-saw mechanism with a variable mass M — M{(j)) 
for the right-handed neutrino. Remarkably the functional dependence M{(j)) has been put in by-hand. It could be 

M{4>) = h<\, or M{<i)) = Moe'^^/f^ with h and / some constants or something else. In the model considered in this 
article, however, we will be able to derive explicitly a form of the neutrino mass as function of the two moduli. In 
our model variable neutrino mass is not put in by-hand. During the expansion of the Universe, the temperature falls 
down till T < rrijj . At that time the neutrinos become non-relativistic particles and an effective potential is generated 
in the form 

Pdark — rTli,Tlu -f- Pdark energy 

(m^). (4.58) 

FNW assumed that pdark is stationary with respect to variations in the neutrino mass. In other words the neutrino 
mass is fixed by a competition between Pdarkenergy and p^. The hypothesis of stationarity implies that 



9 



where 



This last equation determines the temperature dependence of the neutrino mass in a specific dark energy model. 
Adding to this scheme the conservation of energy equation 

p^~iH{p + p), (4.60) 

it is possible to show that 

a; + 1 = = . (4.61) 

Pdark l^i^f^i/ Pdark energy 

Pdark „„\ 

uj = . (4.62) 

Pdark 

Equation (61) inform us that the neutrino contribution to pdark is a small fraction of pdark- Another important result 
can be obtained observing that the continuity equation guarantees pdark ^ i?^'^'^^'"-'. Since n^, ^ R^"^, the neutrino 
mass is nearly inversely proportional to the neutrino density: 

rrii, ~ ~ l/ny. (4.63) 

For further details we refer the reader to references [26] .[271. 

As pointed-out by Peccei [26] , in the FNW scenario ife^ intrinsic scale of dark energy is not obtained naturally since 
the value of the potential in4ife present epoch is fixed by hand: 

V{m1) ^ Timl ^ pDE ~ (2 X lQ-^eV)\ (4.64) 

Thus, the connection of the dark energy scale with the neutrino mass scale is translated into a stabilization problem 
for the field (j). When neutrinos become non-relativistic and m^rii, switches on, is it possible to guarantee the correct 
minimum? 

Remarkably, in our approach we will obtain a neutrino mass dependent on the moduli A and 0. In other words 
the neutrino mass depends on the geometrical configuration of the braneworld. As we already mentioned in the 
introduction, the only dynamic that we contemplate for the scalar fields is a chameleontic behaviour for the radion. 
On cosmological distances, the densities are tiny and the radion can roll on cosmological time scales (the neutrino 
mass is variable with time and possible deviations from a pure cosmological constant contribution can be observed in 
the DE sector); on smaller scales instead, densities are higher and both moduli are simply sitting on the minimum of 
the potential. As a final result, the connection between DE and neutrino mass will be established in the form 

1/1/4 

( U^n ~ 1, (4.65) 

where V is the total moduli potential. 



B. Kaluza-Klein reduction 



Let's consider again the SPINOR action 



SPINOR 



d-x dzV^ i?,- 



m Sgn(Q:) ^'5' 



(4.66) 



We use capital indices A,B,... for objects defined in curved space, and lower-case indices a,b, . . . for objects defined 
in the tangent frame. The matrices 7" = (7*^, «75) provide a four-dimensional representation of the Dirac matrices in 
five-dimensional flat space. The quantity is the inverse vielbein and LUtcA is the spin connection. Because in our 
case the metric is diagonal the connection gives no contribution to the action (|4.66|) . 

The sign change of the mass term under a — > — a is necessary in order to conserve a-parity, as required by the Z2 
orbifold symmetry we impose. 



Using an integration by parts and defining left- and right-handed spinors '^l,r 
written as 



_ 1 



(1 =F 75)^, the action can be 



SspiNOR = d X \J-gA 



1 
2 



(4.67) 



where we considered the redefinition a{z) = e " and we observe that \/— 55 = e \/—gA- The action is even under 
the Z2 orbifold symmetry if ^'^(a;, z) is an odd function of z and ^ir{x^z) is even, or vice versa. To perform the 
Kaluza-Klein decomposition we write 



*L,j?(a;,z) 



(4.68) 
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Because of the Z2 symmetry of the action it is sufficient to restrict the integration over z from ^ to A. Inserting 
the ansatz (|4.68|l into the action and requiring that the result takes the form of the usual Dirac action for massive 
fermions in four dimensions, 

5" = ^ (f x^/-g4,^ipn{x) i^Tpnix) - rrin rp„{x) ipnix)^ , (4.69) 

n •' 

where = tpL + ipji (except for possible chiral modes) and m„ > , we find that the functions fn''^{4>) must obey 
the conditions 

A A 



dz e-/„r (2) f^Piz) = Jdz eVr(^) = S,nn , 

(±9.-m)/„^'«(z) = -m„e'^/„«'^(z). (4.70) 

The boundary conditions fm* i4>) fni'P) — fm* W fnW — 0, which follow since either all left-handed or all right- 
handed functions arc Z2-odd, ensure that the differential operators (±9^ — to) are hermitian and their eigenvalues to„ 
real. Since the equations are real, the functions fn '^{4>) could be chosen real without loss of generality. 

C. The right-handed zero mode: wave function localization 

We now choose boundary conditions such that all left-handed modes are odd under orbifold parity. We specialize 
our investigation to the right-handed massless zero-mode. (|4.70|l now takes the form 

A 

Aze"f^*{z)f^{z) =1, 

(9.+to)/o«(z) = 0. (4.71) 
Let's consider, for example, the limit a 0, then a{z) = e^^^ and we can write 

= 4^—^,^ (4-72) 



where we defined k' = k — 2m and f = f^- 

The zero-mode wave function on the visible brane is thus given by 

jL/ — 2m0 

m = ^F^r^- (4-73) 

Remarkably it is possible to localize the wave function on the hidden brane; in particular for well-separated brane 
we distinguish the cases: 

• k' >0 p{(t>) « 1 k'X » -2m(t) 

• k' <0 p{4>) « 1 ^ e-'"^ « 1. 

D. The Higgs sector and the neutrino mass 

We now describe the realization of a small Dirac neutrino mass in four dimensions, using the localization of the wave 
function mentioned above. Omitting gauge interactions, the action for a Higgs doublet H = (0i,02), a left-handed 
lepton doublet L = [v^, cl) and a right-handed lepton cr localized on the visible brane is 

S = Jd'xV^^^{gZ,,d,Hld.HB-X{\HB\'-vlf} 

+ J d'^x V-.g'"'^^^ {Lsl^df^LB + CRBl^df^eRB ~ {yeLBHBCRB + h.c.) } , (4.74) 

where g[^^^ ^ is the induced metric on the visible brane, and the superscript over the gamma matrices is used to make 
the vielbein implicit. 
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We now introduce a Yukawa coupling of the bulk fermion with the Higgs and lepton fields. With our choice of 
boundary conditions all left-handed Kaluza-Klein modes vanish at the visible brane, so only the right-handed modes 
can couple to the Standard Model fields on the brane. We consider a coupling of the form: 

^-ffind,! [%Lb{x)Hb{x)^r{x, 0) + h.c.} , (4.75) 

where H = i(T2H*, and the Yukawa coupling is naturally of order M~^/^, with M the fundamental Planck s cale o f 
the theory. Rescaling the Standard Model fields and inserting for the bulk fermion the Kaluza-Klein ansatz H4.68|l . 
we find 



x)^^{x)f^{cl>)+h.c.} (4.76) 



where Y5 is naturally of order unity. 

Specializing this formula to the right-handed zero mode we have 



Sy^- J d^x,/^^ii^-j^LB{x)HBix)iPl^{x)f^'{(f>) + h.cj (4.77) 
with a mass for the 4d neutrino that can be expressed as a function of the moduli in the form 

= TO,.(0, A) = \— — , (4.78) 

where fQ^{4>) is given by H4.73|l after the rescaling. 



E. Conformal transformation 

Now we want to study the effect of a conformal transformation to the Einstein frame (g) on the neutrino mass. In 
order to do that, it is necessary to determine the rescaling laws of the various fields in the Yukawa term. We have to 
take into account two different conformal transformations. 

The transformations is completely defined by equations (19) and (20). 

•9^=^9 9tJ.v = Q^gf,^ 

Collecting the two transformations together and remembering [3 — ^"^t"*" we have: 

C"'' = «'(<^)3m^ = fl'(0)g-'.9M- ^ A'~g,,; (4.79) 

5;"''' = a'{X)gf,, = a'{X)Q-'gf.. ^ B'~g^,; (4.80) 

where 

Q^^icoshR)'^^ , S = Q"5F(sinhi?)TTbr_ (4.8I) 

From these formulas, requiring a canonical normalization of the kinetic terms of the fields, we find the following 
conformal rescalings for the fields: 

• CHIRAL ZERO MODE: =^ Q"^/^*^ 



• HIGGS FIELD: Hb =^ He = ^Hb 

• STANDARD MODEL'S SPINOR FIELD: * =^ = (£M)3/2^^ 

Using these last formulas in the Yukawa term, we obtain the "correction" to the neutrino mass induced by the 
conformal transformation to the Einstein frame: 

m.(0,A) = ^^™,(^,A). (4.82) 



12 



V. THE MODULI POTENTIAL 



In the original two-branes model proposed by Randall and Sundrum f|29j.|30j') the brane tensions were tuned in 
order to have Minkowski branes. The separation between the branes becomes a field (radion) in the process of 
compactification to four dimensions. With flat branes the classical potential for the radion is absent: the interbrane 
distance is parameterized by a massless scalar field. Since this field has not been observed experimentally, the authors 
were faced with the so-called "radion stabilization problem". One possible solution is provided by quantum vacuum 
energy of bulk fields: Casimir contribution of bulk fields generates a potential for the radion. The possibility of radion 
stabilization by this mechanism has been considered by several authors using scalar fields ([49], [50], [51j), spinor fields 
f [52. 53] The result is that it is generally not possible to give an acceptable mass toHne radimr and solve the 
hxBiaraiy problem simultaneously. 

When the branes are dS4 or AdSi a classical potential for the radion is present. It is stabilizing for the AdS4 
case and unstable for dS^ branes [38, 54]. The Casimir effect is then a quantum correction to the classical potential 
('[38]). In the dS4 brane case the Casimir contribution has been calculated for conformally coupled scalar fields [55] . 
firinassless spinor fields ([56]) and for massive scalar fields ([79]). The AdSi case has been studied in [57]. 

In this section we will present the classical potential for theTnoduli (Q, R) in the dS^ case. Then we'^oceed with 
the evaluation of the quantum correction induced by a massive spinor field in the bulk. 



A. The classical potential 



In order to write the classical potential it is necessary to review the main characteristics of the supersymmetric 
Randall-Sundrum model (158]) with detuned brane tensions. The detuning process will render our model more realistic: 
since we want to describe^il Universe with a positive cosmological constant, a de Sitter geometry for branes is the 
correct one and detuning is necessary. 

Denoting Ti,2 and Tf respectively the tension of the first (second) brane and the fine-tuned tension, supersymmetry 
imposes the bound |Ti,2| ^ Tj. The detuning process converts Minkowski branes in curved branes. The AdSi solution 
arises when |ri_2| < Tf; the dS4 solution corresponds to |Ti_2| > Tf. 

When one detunes the brane tension Ub TUb, the moduli pick up a potential [34, ^39J : 

^^6(T-l)fc^„,,^ (5.83) 

Taking into account both branes, the classical potential expressed as a function of = InQ and R in the Einstein 
frame is 

Vciass{S,R)^V,ff + W,ff, (5.84) 

where 

V^ffiS, R) = 6(r-l)fc ^.,,,.g/(,^,,.) (coshi?)(^-^"')/(i+2"') (5.85) 

and 

W,ff{S,R) = 6(r-l)fe ^-i2.-s/(i+2.^) (3inhi?)(4-4"^)/(i+2a^) . (5.86) 

Note that T = 1 corresponds to the BPS case. 

As mentioned in [38] the potential is unstable for dSi and we proceed with the analysis of the Casimir energy as a 
stabilizing contribution for the system. 



B. Casimir contribution: the 5-dimensional bag 



We will calculate the Casimir contribution of the bulk spinor on a Euclideanised form of the metric. On the 
Euclidean section the de Sitter branes become concentric four spheres [59], 



ds^ =dz^ +a^{z)dnl, (5.87) 
where dfl^ is the 4-sphere metric. The metric is conformal to a cylinder / x S** [59. 60] . Thus, 



ds^ ^ a'^{z)idy^ +dnl) a{z) ^ {1 ~ Aka^ z)T^ , (5.88) 

where the coordinates are dimensionful and we defined dy — The dimensional length L is given by 

^) ^ t^)^ " - (1 - 4fca2</.)i-ii^]. (5.89) 
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Let's consider a second transformation given by 

y = -Rlnr, (5.90) 
with < r < 1 and i? is a constant introduced for dimensional reasons. In this way H5.88|l becomes 

where dS^ = R^^dil^ is the metric for the 4-sphere of radius one, = a^/r^ is the total conformal factor and r' = Rr 
is the dimensionful radial coordinate (0 < r' < i?). The metric (|5.88|l is thus conformally related to a 5-dimensional 
generalized cone endowed with 4-sphere as a base: a 5-ball of radius R. 

In order to keep track of the geometrical " reconfiguration" of the system under the total conformal transormation, 
let's consider the following formulas: 

• z < z* = where z* corresponds to the position of the bulk singularity, 
' y = —Rlnr. 

The first one corresponds to the assumption that the branes can't fall into the bulk singularity. Since we are interested 
in the description of a low-redshift Universe, we will consider well-separated branes: we assume that the hidden brane 
is close to the bulk singularity {zhidden z* , remarkably this is a general behaviour in the low-energy regime of 
BPS-braneworlds [761 ) while the visible one is located far away in the bulk [z^s « Zhidden)- The remaining formulas 
modify the branesSrlordinate in the following way: 

z ^ -oo y = ^=> r = 1, VISIBLE (5.92) 
z -> z* ^ y ^ +CX) r = 0, SINGULARITY. (5.93) 

Remarkably the total conformal transformation translated the cosmological two-branes set up into a 5-ball. The 
ball radius is connected to the moduli A and (j) by 

^ = -Ine, (5.94) 

where e is a dimensionless parameter that is small for well-separated branes. In this geometrical configuration the 
evaluation of the effective action will be easier, as we now deduce. 

We will start recalling the SPINOR action and studying the effect of the conformal transformation on the spinor 
field \['. The action is 

S = Jd^xJdz^l^E^ '-^j-(dA^^)^+'^^{r,a''^}^ -msgn(«)§vl/| . (5.95) 

In the transformation gAB = 0^9'ab action is rewritten as 

S = Jd^x 1^;^ ' ^ - ^)*' + ^ *'{7°, fi'"}*' - m' sgn(a) *'*'| . (5.96) 

where the primed quantities are referred to the ball metric and they are given by*' = l3'^-^,m' = l3m and E^' = (3E^. 
We will now use a more compact notation and, omitting the prime, we write H5.96|l as 

S = i jd^x (5.97) 

where D = /V + im. 

To proceed further we need to specify boundary conditions for the bulk fermion. Since we are interested in a chiral 
theory on the branes, we choose the "option-L" of [4^, namely: all the left handed modes are odd under orbifold 
parity. In this way the correct boundary condition is — 

z = : P^'iP = 0, z = A : P_7/> = 0, (5.98) 

where P± — ^(1 ±75). To complete the specification of the boundary conditions we remember the existence condition 
for the operator D*: if = 0, this requires that the normal derivative (dy — to)P^_^' = should vanish. In 

summary, we have defined two subspaces in direct sum generated by the operators P±; while on the first space (-) 
we have Dirichlet condition, on the second one (-I-) we impose Robin boundary condition. In other words the bulk 
fermion satisfies mixed boundary conditions (for further details see [52^, ^61, J^ ^' 



14 



In this way our cosmological set up is the 5-dimensional "extension" of the MIT bag model [62]. In these systems 
quarks and gluons are free inside the bag, but they are unable to cross the boundary. This last condition corresponds 
precisely to the mixed boundary conditions discussed above: a chiral theory on the branes corresponds to the condition 
that no quark current is lost through the boundary. The zeta function for massive fermionic fields inside the bag has 
been considered in [6^. Analogous calculations have been developed for massive scalar field ( [64. 65 . 661 ) . Functional 
determinants were otsfcussed in [67. 68, 69, 70, 71, 72, 731. For moduli stabilization with zetaTumTrioTrTegularization 
see also £4, 7^. bf df dfllD' cd 

In the nexTsection we will discuss the zeta function of the system following the approach of [631 . 



C. Casimir contribution: the function 

We start recalling the general set-up for the MIT bag model in three dimensions. The setting we consider first is 
the Dirac spinor inside a spherically symmetric bag confined to it by the appropriate boundary conditions. Thus, we 
must solve the equation: 

if 0„ (r) (r), (5.99) 

H being the Hamiltonian, with the boundary conditions 



1 + i ( ^7 



Cn \r=B. 



=R = 0. (5.100) 



The separation to be carried out in the eigenvalue equation H5.99|l is rather standard and will not be given here 
in detail. Let J be the total angular momentum operator and K the spin projection operator. Then there exists 
a simultaneous set of eigenvectors of H, , J3, K and the parity P. The eigenfunctions for positive eigenvalues 
K = j + 1/2 of is: read 



(5.101) 



whereas, for k = — (j + 1/2), one finds 



i Jj{u}r) Vljim [j) 




(5.102) 



Here lo — \J — m? , A is a normalization constant and fljim{r/r) are the well known spinor harmonics. In order to 
obtain eigenfunctions of the parity operator we must set Z' = Z — 1 in (|5.101|l and Z' = / + 1 in (|5.102|l . In both cases, 
j = 1/2, 3/2, 00, and the eigenvalues are degenerate in m = — j, +j. 

Imposing the boundary conditions (|5.100|l on the solutions H5.1U1|) and (|5.102|1 . respectively, one easily finds the 
corresponding implicit eigenvalue equation. For k > 0, it reads 



E — m 
and for k < 0, on its turn. 



^ ^"^Jj+i{ujR) + Jj{LoR) = 0, (5.103) 



J,(..i?)-^/|^J,+i(c.i?) = 0. (5.104) 



We define the zeta function of the system as 



k 

The power of the method lies in the well defined prescriptions and procedures that we have at our hand to analytically 
continue the series to the rest of the complex s-plane, even when the spectrum Ek is not known explicitly. In particular 
we are interested in the calculation of the Casimir energy for massive spinors in the bag and we will employ zeta 
function to regularize this ground state energy. Briefly we consider 

k 

= -ic(s-l/2)M'^ (5.106) 
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and later analytically continue to the value s = in the complex plane. Here sq is the abscissa of convergence of the 
series, /i the usual mass parameter. 

Let's apply this QCD analysis to the cosmological branes system. When we pass from three to five dimensions we 
have to take into account new degeneracy factors for the eigenvalues, while the implicit eigenvalue equation remains 
formally untouched. In particular the degeneracy d{j) for a spinor field on the ball is 



n + D-2 
D-2 



(5.107) 



where dg is the spinor dimension, D is the manifold dimension (i.e. the dimension of the generalized cone), j = — 1 
and n = 0, 1, 2, .... Since a closed analytical form for the eigenvalues is not available for this system, we will exploit 
the residue theorem (for a pedagogical description of these techniques see [42]'): 



C(^) 



j=3/2 



^ 1 



2m 



jf{kR) ~ J]+i{kR) + ^Jj{kR)Jj+i{kR) 



(5.108) 



Using the method — ordinarily employed in this situation — of deforming the contour which originally encloses the 
singular points on the real axis, until it covers the imaginary axis, after some manipulations we obtain the following 
equivalent expression for C: 



C(s) 



smTTS 



j=3/2,... 



X — In <! z-2j 
oz 



dz 



nR\ i2( N 



2R 



where we defined z = k/ j. As is usual, we will now split the zeta function into two parts: 



(5.109) 



N 



as)^ZN{s)+ E Ms), 



(5.110) 



namely a regular one, Zn , and a remainder that contains the contributions of the first terms of the Bessel functions 
I^{k) as I/, fc — > cxD with v/k fixed. The regular part of the zeta function is: 
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N 
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(5.111) 



where rj = y/l + z^ + ln[z/(l + vT+z^)] and t = l/vT+z^. After renaming mR — x, the relevant polynomials 
(see appendix A) are given by 
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The number N of terms to be subtracted must be high enough in order to absorb all possible divergent contributions 
into the groundstate energy. In our case N=5. The asymptotic contributions Ai{s), i — —1, 5 are defined as 
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This completes the description of the zeta function. In the next section we will consider the divergent contributions 
in the vacuum energy. After that we will move to the renormalization process. 

D. Casimir contribution: Residue [(^( — 1/2)] 

As a first step towards renormalization, let's start considering the divergent terms in the groundstate energy. By 
construction they are contained in the functions Ai{s). We will discuss them separately following the approach of 

1. Residue A-i{-l/2) 

The first asymptotic contribution reads 

sin(7rs) 



of which we need the analytical continuation to s = —1/2. With the substitution t = {zj /RY — rn? ^ this expression 
results into the following one 

00 °° 

^ sinM ^ I dt^Z^|Vj2 + i?2(t + ,„2)_j| 

TT J t + L J 
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where the Mellin integral representation for the single factors has been used. As we see, the /^-integral is well defined. 
In order to deal separately with the two contributions in the /^-integral, we introduce a regularization parameter 5. 
In this way A_i(s) can be written as 

= Ihn [A\{s,5)+A\{s,S)\ , (5.116) 



with 
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In A^_i{s, 5) two of the integrals can be done, yielding 
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Exploiting the identity 
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it is possible to rewrite A^_^{s, 5) in the form 

oo 

nl-2(5 °° /■ 

A\{s,5) = -^-^r(. + <5-l/2) E d{j) jdxx^-' 
Adding up (|5.117|l and IjS. 119(1 yields 
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where in the last step a partial integration has been performed such that the x-integral is well behaved at s 
This is a form suited for the treatment of the sum, which is performed using 



oo / 



1^=1/2,3/2, 



1 ?Tjr(£_n^-^ 2«+2 

2 r(.) 



+ (~1)"2 / dv 



,2n+l 



l + e2 



+ (-l)"2cos(7rs) / diy 



,2n+l 



l + e2 



-1 -1 

X 



for odd powers of v and 



i'=l/2,3/2,. 



i r(n + l/2)r(.-n-l/2) 

2 r(.) 



-(-l)"2sin(7rs) / rfi^ 



1 + e 



27ri/ 



(9' 



- 1 
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for even powers of v. This formula as well as IjS. 121(1 can only be applied for Re{s) < 1, otherwise the integral diverges 
at the integration limit v = x. 

Since the prefactor in (|5.12U|I contains a pole at s=-l/2 an expansion of the integrand is necessary and remembering 

that d{j) ~ {2/3){j^ + ^ I ^ §)) our final expression for the residue is found to be 

. , , 457 l7m^R m^R^ m^R^ ,^ 

res -1/2) = \ V . 5.123) 

^ ' ' 2419207ri? 28807r 1447r ISOtt ^ ' 



2. Residue Ai{~l/'2:), i> I 
Let's start remembering the form of the asymptotic contributions for i > 1: 



Ms) 



sin(7rs) 



We need to introduce some notation as follows: 



dz 



■R/i 



2i 



d D^{t) 
dz 



a=0 

(where the relevant coefficients Xi^a are explicitly listed in H5.112() and 

!v=l/2,3/2,... \ ^ 



(5.124) 



(5.125) 



-s-b 



With this notation we write 



Ms) 



sin(7rs) d{j) 



2i 



E ^E-- 



TT ^ — ' T 

j=3/2 a=0 

2i 



■R/j 



I ( ^ 1 - 



oz 



a-\-i 



sin(7rg) ^ J x.^^m-^'^ r(s + ^)r(l - s) 

a— I ^ 2 / 



i=3/2 



mR 



sin(7rs)r(l — s) 2 

^TT 3 ^ 



Xi^a r(s + 2 ) 

^ 1 (mi?)"+^ r(2±i) 

a— k ^ 2 ^ 



/(s;3 + a;^) + -/(s;2 + a;^) 



1 , i + a 3 i + a, 

--/(s; 1 + a; ) - -/(s; a; — 

where we performed the z-integration using the identity (\A2]) 



dz 



mR/v 



±r = ^^-2. " r(. + t)r(i-.) ^,^ 
2(mi?)" r (1 + f ) 



mi? 



(5.126) 



(5.127) 



and we exploited the vanishing of the degeneracy factor for j = 1/2. In order to evaluate the /(s; a; b) for the relevant 
values in s = —1/2 we can use the following formulas: 

• the recurrence formula 

/(s; a: b) = (mR)^ [/(s; a - 2; 6 - 1) - f{s; a - 2; b)] ; (5.128) 

• formulas Emu and Eini valid for SRs < 1; 

• an extension of formulas 15.1211 and 15 . 1 2 2| t o dis < k + 1 < n + 2 obtained using partial integrations: 

2\"' n!r(s - n - l)a;2"+2 



E -'"^^ (i + (- 

i/=l/2 



2r(s) 



r(s) 
r(s - /c) 



+2(-l)"cos(7rs) 



r(s) 



7 2 \ ^ 

d X \ V 



,2n+l 



dv2v 1 + e2 



2r(s) 



1- (- 



a; \ V 



2n+l 



dv 2v I 1 + e 



2iri^ 



- 1 



-s+fe 



(5.129) 
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for odd powers of v and 



-s_ r(n+l/2)r(s-n-l/2) 
'^x' ) ~ 2r(s) 



2n+l 



for even powers of v. 



In this way it is possible to construct an explicit form for the asymptotic contributions and, in particular, our result 
for the residues is given by 

_ _ 17m 17 _ m^R nt'B? m^W' m^i?" 

res i(- / ) - ~28807r 115207rJZ ~ 2887r 247r ^ UAtt 127r 

res ^.(-1/2) = -(6^+i^)-^(y4+i4?) 

^512 487r'^^32 47r^ 16 

. / , /I 7,1,1 97 , 2^/ 25 565 , 

res A3(-l/2) = -m(— + ^) + -(^ + ^^^) - + 

3^2/ 3 41 , .3,1 2 . m^R^ 
-- ^^32 + 120^) -"^^(16 + 9;^)-^ 

. / . / 63 11 , 1 , 35 13 , 2r,/ 13 1 , 

re. A,(-l/2) = -m(— + _) + -(_ + _) - + -) 

. / , /^N .61 37 , 1 , 23 31 , 2../ 23 359 , 

res A5(-l/2) = + ^) + ^(y^ - ^^^n^^ +"^^^768 + 3780^^^ 

5. Residue Ao{-l/2) 

This calculation will exploit the basic formulas we showed above and the Mellin integral representation. With the 
substitution t = {xvjaf' — rv? we write 



Ao(a) = ^ f: r J(;±\ 



/ dz 



1 sin(7rs) 
2' 



00 00 

= £ Idt t-^ f da e-«(*+-^) (5.132) 

00 



X 



where in the last step the Mellin representation for the single factors has been used and, in particular, we exploited 



V f + R'it + -^')- l =Jfdp p-^/^ r _ . (5.133) 

In this way wc can closely follow the calculation discussed previously for Res 1/2). The asymptotic contri- 

bution is rewritten as 



= -2^T(^^.g/W^|[,-2 + (^^)2]«+V2 
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. + imRrj\. jjt:^^^^,^} (5.134) 



and the relevant residue is given by 



E. Casimir contribution: renormalization 

For the discussion of the renormahzation let us look for the divergent terms in the groundstate energy. Adding up 
the contributions calculated previously, the residue for the zeta function is 

_ 4m 41 2m'^R 23m^R'^ 7m^R^ m^R* m^R^ 

Res C(-l/2) - ^ + 10395^^ " ^ ^ + + (5-136) 

These terms form the minimal set of counterterms necessary in order to renormalize our theory. 
We are led into a physical system consisting of two parts: 

1. A classical system consisting of a spherical surface ('bag') with radius R. Its energy reads: 

Eaiass = aR^ + bR'^ + cR^ + dR^ + eR + f + -^. (5.137) 

R 

This energy is determined by the parameters: a, 5, c, d, e, /, g. 

2. A spinor quantum field 'i'{x) obeying the Dirac equation and the MIT boundary conditions H5.100|) on the 
surface. The quantum field has a ground state energy given by Eq, Eq. (|5. 106(1 . 

Thus, the complete energy of the physical system is 

E^Eciass+Eo (5.138) 

and in this context the renormalization can be achieved by shifting the parameters in Edass by an amount which 
cancels the divergent contributions. 

First we perform a kind of minimal subtraction, where only the divergent contribution is eliminated, 

1 ml' 1 

0^0 



3607r s 1207r s 

7m^ 1 . . 23to3 1 



d^d--—- (5.139) 



ISOtt s 6307r s 

m? 1 2m 1 

/ f + 

457r s 3157r s 

41 1 



9^9 



207907r s 



The quantities a — {a, 6, c, d, e, /, g} are a set of free parameters of the theory to be determined experimentally. In 
principle we are free to perform finite renormalizations at our choice of all the parameters. 

We find natural to perform two further renormalizations. First it is possible to determine the asymptotic behavior 
of the Ai for to — > cx). The finite pieces not vanishing in the limit m — > oo are all of the same type appearing in 
the classical energy. Our first finite renormalization is such that those pieces are cancelled. As a result, only the 
"quantum contributions" are finally included, because, physically, a quantum field of infinite mass is not expected to 

fluctuate. The resulting Ai will be called A^[^^\ 

Concerning Z we have constructed a numerical fit of Z by a polynomial of the form 

6 

-P(m) = ^qto", 

1=0 

and then subtracted this polynomial from Z . As explained above, this is nothing else than an ulterior finite renor- 
malization. The result will be denoted by Z^'^'^"-\ 
Summing up, we can write the complete energy as 

£^ = i^c/ass + (5-140) 

where Edass is defined as in H5.137|l with the rcnormalizcd parameters a and E'q'^'^"-' = -|- Y^\=^i A^['^"'\ 

Figure ^ shows the quantum contribution e'^'^"'^ as a function of the bag radius. Remarkably the Casimir con- 
tribution of the massive spinor is stabilizing: the energy exhibits a minimum corresponding to a stable geometrical 
configuration. 

In the next section we will discuss the parameters of the model and the stabilization mechanism in greater detail. 
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VI. MODULI STABILIZATION AND PARAMETERS SPACE 



Is it possible to achieve a realistic model with a "natural" choice of the parameters? To answer this question it 
is necessary to remember that our purpose is to stabilize the entire geometrical configuration of the system which 
amounts to stabilizing both moduli. As we already said, it is not possible to introduce in a model ultralight scalar fields 
with a generic coupling with matter, since phenomenological constraints must be faced. In this model the two moduli 
interact with matter in a known way. In fact, we could parametrize the strength of this coupling, since we know 
explicitly the conformal factors (4.78), but our approach will be different: dangerous couplings are kept under control 
because we are giving a (large enough) mass to the moduli. It seems now worthwhile to focus on the stabilization 
mechanism. 



A. R-modulus stabilization: the chameleon mechanism 



As we already mentioned in section 3, the possible interpretation of the R-modulus as a chameleon field has been 
investigated in [48] . Following their approach wc will add a run-away bare potential for the radion in the form: 



Krf<i(i?) = Af*eV">'"y (6.141) 

where M = 10~^ eV in order to trigger acceleration now. It is not a quintessence potential since it does not converge 
to zero at infinity. It should be considered as a very flat potential which may appear from some non-perturbative 
physics for the radion. 

For the potential above, the minimum of the effective potential 

VeffiR) - VaddiR) + Pm(l + ^R^) (6.142) 



is given by 



i?"+2 = n i —\ ^'^ddjRmzn) ,g 



(see formulas 13. 5514. 81|) . 

For Vadd = O(M^) this leads to a cosmological value for the radion field given by 



R'^'^Oi — ] . (6.144) 



Notice that Roo = 0(10^^*^) for n = 1 implying that solar tests are automatically satisfied. 
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The mass of the radion at the minimum reads 

M _ 1 



n{n + 2) + n\ — r — 



(6.145) 



The interesting regime is obtained for small distances, then the matter density is high and m H, the field "sits" 
on the minimum of the effective potential (i.e. the minimum acts as an attractor): the radion is stabilized. 

Although the chameleon mechanism does work also without runaway potentials [771 , the presence of the potential 
I6.141l is crucial in our model. Had we considered only the classical potential 15. 851 foFrhe radion, we would have been 
forced to face the following problems: 



It is not possible to achieve the competition between the potential {Vdass ~ Cosh{R) i+2q^ ) and the matter 
coupling (see formula [4. 8 



• The value i? = for the radion is not forbidden by the classical potential. The falling of the hidden brane into 
the naked singularity (i? = 0) must be avoided if our intention is to maintain a supergravity approximation [351 . 



It seems noteworthy that the Q-modulus is not a chameleon since it is not possible to achieve the desired competition 

between the potential (VdassiQ) ^ Q^^) and the coupling to matter {A{Q) ^ Q~^). To stabilize the Q-modulus 
we suggest to exploit the Casimir energy of the bulk spinor evaluated in the previous section. 



B. Q-modulus stabilization: the Casimir energy 



The stabilization of the system is based on the presence of (A) two scalar degrees of freedom (the positions of the 
two branes are parametrized independently) and of (B) a cosmological attractor that pushes the hidden brane towards 
the singularity. 

As mentioned above, the Casimir energy of the massive spinor guarantees a stable geometrical configuration for the 
5D bag. The connection between the bag and the braneworld scenario is the conformal transformation (|5.9U|) . Had 
we had only one scalar field in the low energy description of the model, than only the distance between the two branes 
would have been physically relevant (consider the radion in Randall-Sundrum two branes model) or, equivalently, 
it would have been possible to shift the zero of the fifth dimension: moving the two branes away in the bulk while 
keeping fixed the distance between them. In this case the conformal transformation H5.9()(l would not have mapped 
the two branes set up into a ball, since only the bulk singularity corresponds to the centre of the ball. 

The second crucial ingredient is the presence of a cosmological attractor that pushes the hidden brane towards the 
bulk singularity [761 . The visible brane evolution becomes independent of the hidden brane and the braneworld system 
is described by aT7n(?/e-scalar-tensor theory. For this reason we will assume a very small value for the R-modulus (the 
hidden brane is close to the singularity/bag centre) in a late-time cosmology we are interested in. In this way we can 
" forget" about the R modulus (it is fixed at small values and stabilized by the chameleon mechanism) if the efficiency 
of the attractor is high enough. The Q-modulus corresponds to the remaining scalar degree of freedom (position of 
the visible brane/bag surface) and the last step in the stabilization process will be to add the Casimir energy as a 
stabilizing contribution for the bag surface. 

To proceed further, it would be necessary to study in full detail the role of the conformal transformation to the 
Einstein frame (i.e. the cocycle function, for a pedagogical treatment see [42]) and write the full moduli potential. 
However, we will consider a reasonable value for the Q-modulus {Q ^ I) for ihe discussion of the neutrino mass. 



C. Parameters Space: Neutrino Dark Energy 

Here is a possible choice of parameters. We will proceed stepwise: 

1) set the fundamental scale of gravity as 

kl = M/w, (6.146) 

IQ^^GeV. (6.147) 

2) Choose ^ 10~^ in order to face with phenomenological constraints ([37, 761). 

3) Set the tension scale and the 5D neutrino mass near the fundamental seme or gravity: 

/fc - m - lO^^GeV. (6.148) 

4) The Newton's constant as measured by Cavendish experiments in the visible brane is given by (see also |76|): 

^ {lO^^'GeVf. (6.149) 



A:/c|(l + 2a2)a2(0) g^Q 
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For the screening potential of the form H6.141|l the mass scale has to be chosen such that M « 10~^ eV if the radion 
must be a dark energy ca ndida te. This tuning is no more than that required by a cosmological constant. 
The classical potential (|5.85|l shows a scale given by 

V^iass ^ ^ ' 6(T - l)Ml,,anck- (6.150) 

The correct dark energy scale is recovered imposing the usual (cosmological constant) fine-tuning T — 1 ~ 10^ 

On cosmological distances, if we consider rt = 15 in the potential we obtain i?oo ^ 10^^''. In this way the neutrino 
mass is strongly suppressed by the large separation of the branes. In fact, choosing 



k 
20 



(6.151) 



and evaluating the neutrino mass in the Einstein frame for Q ^ 1 and i? ~ 10 (corresponding to the minimum of 
the moduli potential after the discussed stabilization), we can write 



= Y5f^{q})vi,^^^ ^ 10-^eV. (6.152) 

Remarkably, it is possible to obtain a meV-mass for neutrinos with a reasonable branes configuration in five 
dimensions. Happily, in the minimum (Q, R) of the total potential we can write 

'-i^'S, = fiT, h, k, m, a, M, n) ~ 1, (6.153) 

without fine-tuning the parameters (with the only exception of the usual cosmological constant fine-tuning, a proper 
choice of n = 15 in the potential and a small value of a). With this last configuration the radion mass on cosmological 
distances is m ^ lOiJ, the field can roll and the neutrino mass is variable with time. Also, notice that the neutrino 
mass is space-dependent: this can be a source of phenomenological consequences. 



VII. CONCLUSIONS 



In this paper a braneworld model for neutrino Dark Energy has been presented. The model is a supersymmetric 
extension of the Randall-Sundrum two branes model. After supersymmetry breaking the system becomes dynamical 
and a (classical) moduli potential is generated. 

In this article we have achieved the following results: 
1) In our model a 4D neutrino mass as a function of two scalar degrees of freedom (moduli) has been explicitly 
calculated rather than inserted by-hand. 

2) The moduli stabilization issue is addressed with a de Sitter geometry for the branes. The chameleon mechanism 
is exploited to stabilize the radion following the approach of [48] . As far as the stabilization of the remaining modulus 
is concerned, we suggest to exploit the Casimir energy of massive neutrinos as a stabilizing contribution. 

3) A direct relationship has been established between neutrino mass scale and Dark Energy scale. We stress 
again that the only dynamical behaviour associated to moduli is encoded into the chameleon field R. Namely, on 
cosmological distances, the densities are tiny and the radion can roll on cosmological time scales, while on smaller 
scales, the chameleon acquires a large mass (radion stabilization). The neutrino mass is variable with time on 
cosmological distances and possible deviations from a pure cosmological constant behaviour can be observed in the 
DE sector (we will further discuss these issues in a future work) . The neutrino mass is also space dependent and this 
can imply phenomenological consequences. 

It seems to us that this model suggests a stronger connection between DE and neutrino physics. To the best of our 
knowledge it is the first attempt of connecting neutrino mass to Dark Energy from the standpoint of brane models. 
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APPENDIX A: D-POLYNOMIALS 

In this appendix we collect the basic formulas necessary to obtain the explicit form of the polynomials 15.1121 We 
start remembering the uniform asymptotic expansion of the Bessel functions Iu{k) as i^, fc — > cx) with v /k fixed. One 
has 
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1 e"'' 

(1 + 



fe=l 



with t ~ + and 77 = \/l + + ln[z/(l + \/l + z^)]. The coefEcients are determined by the recurrence formula 



Uk+i{t) = ^t^l - t')uUt) + i / dr (1 - 5r2)ufe(T), 



starting with uo(i) = 1. As is clear, all the Uk(t) are polynomials in t. 

Since boundary conditions involved Bessel functions and their derivatives, we need also the expansion 



/27W z 



k=l 



with the ffc(i) determined by 

vk{t)^uk{t) + t{e -\) 

Exploiting the previous formulas the relevant asymptotics can be found. With the notation 



Si 



1+E 



Uk(t) 



k=l 



1+E 



Vk{t) 



k=l 



, mR = X 



it reads 



ln[/|(zj)(l + 1 - ^) + + - j)I,izj)mzj)] 



re2J"(l + z2)i/2 1 2x, z2 „ , 2, , 1 ,1 

1^ 27rjz'^ z^j 1 + z'^ ] (1 + z'^)-^/^ J 



In \ C V 2(1 - t))^ ' — 



27rjz2 2(1 -f) 

In this way we are led to the definition of D-polynomials as follows 



E 

fc=i 



= In I ^^[(1 - + + 2^^E,S,] 



2(1 - i) J 
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